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Abstract 

The problem of phase retrieval is a difficult one which remains far from solved. 
Two homometric sets are always connected by way of a convolution product by 
some spectral unit, though not necessarily in a unique way. Here we elucidate 
one small aspect, the subgroup of spectral units with finite order. Its elements are 
completely characterized by relations between their eigenvalues. This sheds some 
light on the beltway problem. 
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1 Introduction 

In this paper we fix some integer n and consider vectors in <C n as maps from 
the cyclic group Z n to C: a = (a(0), . . . a(n — 1)) = (a ■ ■ ■ a n -i)- Usually, 
such a vector can be the characteristic function of a subset of Z 
put a = (1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0) for subset {0, 4, 7}. Two vectors a, b are 
homometric if equivalently 

• They have the same interval distributions (i.e., histogram). 

• Their Fourier coefficients have the same magnitudes. 

This problem originated with diffraction patterns, as illumination of such patterns 
is essentially proportional to the amplitude of a Fourier transform. If we know 
that a is homometric to some given set b then finding a exactly is known as the 
phase retrieval problem, because already the amplitude of a's Fourier coefficients 
is known. The problem is that a is not necessarily congruent to b, as seen from the 
example: 
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Example 1 Consider the characteristic functions of {0, 1, 4, 6} and {0, 1, 3, 7}, 
subsets ofZ 12 i.e. a = (1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0), b = (1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0). 
The sets are not congruent but are homometric: for instance the inter\'al distribu- 
tions are the same, e.g. (4,1,1,1,1,1,2,1,1,1,1,1). Alternatively the Fourier coeffi- 
cients have magnitudes (4, V2, 2, y/2, 2, y/2, 2, a/2, 2, V2, 2, y/2). 

In a paper with William Sethares (T| about musical scales or rhythms expressed 
in terms one of another, we introduced scale matrices, which are circulating ma- 
trices with a given first column. 

Definition 1 AI £ M n (C) is a circulating matrix iff its general coefficient is 
rriij — Q-i-j (where i — j is computed modulo n). 

Hence any circulating matrix AI can be written M : aoI n +Qi J+ ■ ■ ■ a n -i J" -1 

/ o o ... 1 N 

' 1 0... 

where the vector a = (ao, Qi, . . . a n -i) is given and J = 

\ 6 ... io/ 

The scale matrix associated with the scale vector a — (ao, oi, . . . a n -i) is 
precisely 



/ ao On-l 

ai ao 



M(a) = a I n + aiJ ■ 



jn—l 

. a n -\j — 



CLn-1 



ai 



ai ao J 



\a„-i 

All these matrices form a commutative algebra C„ . 

The matrix product M(a) x M(b) corresponds with the convolution product 
of a and 6, as seen from the first column of the product matrix: 

71—1 TO — 1 

M(a) x M(b) =M{a-k b),a*b= a k b- k , ^ a k hi-k, ■■■) 

k=0 k = 

where all indexes are computed modulo n. 

The connection with Fourier coefficients is straightforward: 



Proposition 1 We define the Fourier matrix Q, by its (j, k) th term, equal to ■ 

Notice that 'fi = f2 and fl^ 1 — Q. Then for any a G C n ,Q~ AI(a)Q is equal to 
the diagonal matrix F a with coefficients 



a 2ijkif/n 



F a (0) = J2a k) ... F a (j) = «fee 2ljW ", . . . F a {n-1) = ^ a k e 2i(n - 1)kn/n = J2 

fc=0 fc=0 fc=0 fc=0 

i.e. the Fourier coefficients of a. 

We retrieve from there the famous relation T(a ★ b) = J-(a) x In more 

abstract terms, three algebras are isomorphic: 



-2ikir /r. 



1 This explains why the equality of the interval contents is equivalent to the equality of amplitudes of 
Fourier coefficients, as the interval content is computed from a * (—6). 



3 



1. the algebra of columns C n , identified with the maps from Z n to C, and with 
the convolution product; 

2. the algebra of circulating matrices, and 

3. the algebra of diagonal matrices. 

Homometry is then readable on the matrix with the magnitudes of the Fourier 
coefficients, which is (squared) 

n~ 1 M(a)n. t n- 1 M(a)Q. = fi _1 M(a)*M(a)fi 

For real valued vectors a we have of course M(a) = M(a). In other words, a and 
b are homometric iff a * a* = b * b* i.e. M{afM(a) = M (fo)'M(6). (NB: a* is 
derived from a by retrogradation and conjugation) 

Another way to look at it is to consider a diagonal matrix D such that 

F b = Q _1 M(6)n = D!i _1 M(a)!i = DF a 

where the diagonal elements of D are of amplitude 1 for a, b to be homometric. 
We get in return a characterization: 

Proposition 2 a, b are homometric if there exists a matrix U such that 

• *U = IJ- 1 

• U€C n 

• B = UA. 

Equivalently, we have b = u * a with u * u* = (1,0, 0,...) (the more usual 
definition, see J2|). Such elements satisfying u * u* = (1, 0, ... ) = <5o, i.e. 
homometric to So, are called spectral units. 

The difficult issue there is: where do the coefficients of u (or U) lie ? Ab- 
stractly, the group of all possible C/'s is a torus, as each eigenvalue can be any 
complex number on the unit circle (the eigenvectors are always the same, since 
they are shared by all circulating matrixes). In practice (for instance in music the- 
ory), we look for integer, or at worst rational, coefficients. Rosenblatt proved (his 
Thm. 3.1, (2)) that 

Proposition 3 For any given subfield K ofC, closed under conjugation, for a,b £ 
K n , a and b are homometric iff there exists a spectral unit u £ K n satisfying 
b = u-k a. 

For instance, u — (0, 1, 0, ... 0) gives rise to a permutation circular matrix U. 
As a spectral unit it realizes the translation of a subset by 1. It is worthy of note that 
symmetries of the cyclic group Z„, mapping the characteristic function of some 
subset S to the characteristic function of, say, — S, are much more complicated in 
terms of spectral units. For instance the spectral unit connecting 

1(0,3,7} = (1,0,0,1,0,0,0,1,0,0,0,0) 
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and 

l7-{0,3,7} = 1(0,4,7} = (1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0) 

(minor and major triads in music theory), is 

(7,4,-2,l,7,4,-2,l,-8,4,-2,l)/15. 

In the context of this paper, it means that B = UA where (taking henceforth 

K = Q) 

• U is unitary: '17 = U" 1 , 

• Us C n , and 

• U is rational valued: U £ M n (Q) 

NB: the result is easy if A (and hence B) are non singular, because U = BA~ l . 
We call such matrices spectral units again. 

Unfortunately, as the example above with the major and minor triad shows, 
some spectral units are not simple - the one in the example has infinite order. This 
paper enumerates all (rational) spectral units with finite order. Since even in the 
unit circle considered as a group, both elements with a finite order and rational 
elements form infinite subgroups, it is notable that these spectral units are finitely 
many. To be specific we will show three theorems. Remember we are dealing with 
rational spectral units with finite order only. 

Theorem 1 If U is a spectral unit with finite order and n is even, then all its 
eigenvalues are n th roots of unity. Ifn is odd, then the eigenvalues are either n th 
roots of unity or their opposites (i.e. 2n* root of unity). 

This stems from a more precise condition 

Theorem 2 IfUis a spectral unit with finite order and n is even, then for all k 
coprime with n and any Fourier coefficient ( = eigenvalue of U )£j ,j / 0, one has 
£kj = For j = we have £o = ±1. 

This condition stands even when n is odd, except when £j is a e ( 2 Pj + 1 ' l7r / n 
and k even, when we have £jy = — £j . 

For instance for k = — 1 this gives the condition that the last Fourier coefficients 
must be the conjugates of the first ones (thus ensuring that U is real valued). 

These conditions enable to specify exactly what are the possible values of the 
Fourier coefficients £, and hence to retrieve all possible spectral units with finite 
order (by inverse Fourier transform). We will see some examples in the next sec- 
tion. A general result can be stated which gives the precise cardinality of the group 
of such spectral units, as it gives directly the decomposition of this (abelian) group 
into a product of cyclic groups. 

Theorem 3 Any spectral unit with finite order is defined by the values of the subset 
| n} of its eigenvalues. The possibilities are listed infra: 

• Co = ±i; 
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• When n is odd, for all j \ n, OR —£j is any power ofe 21 ^'". 

• When n is even, £j is any power of e 2 ^ 71 ^ ifn/j is even, or any power of 

Unfortunately, many spectral units are not of finite order (even in simple cases), 
like (0, 1/2, 0, 1/2, 1/2, -1/2, 0). Still the results above may perhaps enable to 
compute all spectral units with, say, small denominators, which occur in practice 
for homometric subsets of Z„. 



2 Examples 

Example 2 Let us elucidate the group of units when n = 12. Let u be a spectral 
unit with finite order, and £o, • • • £n its Fourier coefficients, i.e. eigenvalues of the 
associated matrix U. From theorem 2 above, the relation = £jk is satisfied for 
all four values of k = 1, 5, 7, 11. 

• There are no conditions on £i which is any 12 th root of unity; its value 
specifies £5 = £ 5 and similarly £7, £11. 

• £2 must be a power of e 2x2 ' l7T / 12 — e 17 '^: without making use of Thm. 3, 
we can see that if we set £2 = e 2lrra ^ 12 where a 6 Z12, then £f = £10 = 
£_2 = £2 > hence 5a = —a mod 12, i.e. a is even. This determines also 
£10 = ? a - 

• Similarly £3 is a power ofi = e 1 ^ 1 " 2 : if we set again £3 = e 2 " l7Ta / 12 we g e t 

£3 = £9 = £-3 = C^ 1 i.e. 7a = —a mod 12 

meaning that 8a — i.e. 3 divides a. We have £9 = £ 3 . 

• ?4 gi ves - fne slightly more complicated case: it must be a power ofe ln ^ 3 , just 

T.i t j . f 2i7r/3 2x4i7r/12 ■ r <- 2iwa/12 

like £2, and not a power of e ' = e ' : ;/ £4 = e 

£4 = £g = = £4 1 i.e. 5a = —a mod 12 

Here also, we find that £s = £4- 

• fo ± 1 and £e = £_6 = £g 1 is a 12/6" 1 roo« o/7, ;'.e. f;e = ±1. 

7b conclude: £1 is any 12" 1 root of unity, while £2, £3, £4 are limited to subgroups, 
£0 and ^6 = ±1- 77ie structure of the group is then TL\% x (Ze) 2 x Z4 x (Z2) 2 , 
wifn 6,972 elements. 

Example 3 7afe n = 7. Then apart from £0 = ±1. all Fourier coefficients are 
powers of '£1: 

• 7f £j = 1 /nen £2 = £1 and more generally = , fe = 2 . . . 6. 



6 



• If £i 7^ 1 — Ci 4 then we play the same thing with — £1.' £2 = — £1,6 = 

a /^6 = r x (= 6 3 ),6 = 6 _1 = &6 = c 1 - 

In the exponential notation, we get more simply 

(2fc + l)i7r 2(2fc + l)i7r 6(2fc + l)i7T 

(6,6,6,6,6,6) = (~e 7 ,-e 7 , e ? ) 

When 6 = +1 gives spectral units of the form 

u = (2/7, -5/7, 2/7, 2/7, 2/7, 2/7, 2/7) 

or a circular permutation thereof; if£o = — 1 ifen we obtain opposites of 
circular permutation matrices, e.g. u = (0, —1, 0, 0, 0, 0, 0). These spectral 
units are exactly the opposite of the ones obtained in the first case, when 

a = 1. 

The whole group has only 28 elements (its quotient by ±1 is the cyclic group of 
14 th roots of unity). 

Example 4 We end these examples with an odd composite number: n — 15. As 
usual fo = il. 

• Coefficients k € Z*5, are deduced from £1 = £: 

- //£ 1S = 1 (e.g. £ = e 8i?r/15 J fto?n 

£2 = £ 2 ,6 = £ 4 ,6 = .f,6 = f 8 , 

€ 11 =f 11 ,6s = f ls ,Ci4 = f 14 = r 1 =f 

e.g. Cfc = e 8fci,r/15 /orA: = 2,4, 7,8,11,13, 14. 

- /f^ 1 = — 1 instead, then the list is the opposite of the one found in the 
first case, e.g. £ k = -e skiw/15 for k = 2, 4, 7, 8, 11, 13, 14. 

• Similarly, £g, 6, 62 ^W" 6 powers 0/^3 for fneir opposites), which must be a 
5 th root of 1 (or -1): for example, 

/- 4iir/5 /- 8i7r/5 /- 16i7r/5 6i7r/5 ^ ^8 32i7r/5 2i7r/5 

6 = e , 6 = e ' , €12 = e = e , £ 9 = £ 3 = e = e 
or (when £3 is a fifth root of -1) the opposite: 

6 = — e , 6 = -e ,Ci2 = -e ,6 = -e 

• 7"we same goes for £5 = -|- e 2lfc,, 7 3 an< ^ £ 1Q _ 

More concisely: £1 is any 30 th roof 0/ unity, £3 any 10 t ' 1 roo«, £5 any 6 th root. 
The group has 2 x 30 x 10 x 6 = 3, 600 elements. 

Example 5 Say a, b are homometric. If matrices A, B are non singular, then there 
is one and only one spectral unit u, with matrix U = BA~ X such that B = UA 
(or equivalently b = a-ku). But sometimes eigenvalues of A, B are nil. ^2ty chooses 
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to give 1 as eigenvalues for the corresponding eigenvectors, but other choices are 
possible, if we respect the conditions found in the theorems above. 

An example issued from music theory is a — (1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1), 6 = 
(1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1), two melodic minor scales. The Fourier coefficients 
with indexes 2 and 10 are nil. 

• Using Rosenblatt's choice, we take £2 = ^ti(2) = •T-'u(lO) = £10 = 1 
(the other Fourier coefficients are determined by J- U (k) = J~b(k) / J- a (k)). 
This yields u = (0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0). Musically this means that A 
minor is transposed from C minor by a minor third. 

• We know from Thm. \2\that £2 is some power ofe™^ 3 , £10 being its conjugate 
or inverse. This yields five other possible units, e.g. 

u= 1(1, 0,-1, -1,0, 1,1, 0,-1, 3,0,1) or 

u = ^(1, 2, 1, -1, -2, -1, 1, 2, 1, 11, -2,-1) or 

u = 1(2, 1, -1, -2, -1, 1, 2, 1, -1, 4, -1, 1) or 
6 

u = -5(1, -1, -2, -1, 1, 2, 1, -1, -2, 11, 1, 2) or 

u= 1(1, 1,0,-1,-1,0, 1,1,0,3,-1,0). 

In a way this can be interpreted as other, hidden symmetries between those two 
musical scales. 

3 Proofs 

3.1 Proof of Thm. U 

Throughout, U is a circulating matrix which is unitary (U^ 1 = *[/), is of finite 
order: U m — I n for some m, and has rational elements. Hence its eigenvalues 
are of magnitude 1 (they are m th roots of unity), and as discussed above U diago- 
nalizes into Diag (£0, £1, ■ • ■ where the eigenvalues £j are also the Fourier 

coefficients of the first column of U, seen as a map from Z n to C. 
We prove an alternative form of Thm. Q] 

Proposition 4 All eigenvalues of U are n th roots of unity for even n, and 2n th 
roots of unity for odd n. 

As we assume U has finite order, all these eigenvalues are roots of unity. Moreover, 
as U is a polynomial U = P(J),P € Q[X] in the matrix J, whose eigenvalues 
are the n th roots of unity, the eigenvalues of U are polynomials in these roots: 
£ fc = p( e 2ik7v/n ) , hence lie in the field Q[e 2in/n ] = Q„, often called cyclotomic 
field. We need the following 
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Lemma 1 Let £ be a ml h root of unity belonging to the cyclotomic field Q n . 

(£™ = 1 when n is even, 
£ 2n = 1 vv/jerc n is odd. 
In other words, if Q m C Q n then m is at most n or 2n, according to whether n is 

even or odd| 

Let £ be such a number (like any eigenvalue of U). Let m be the order of £, 
i.e. the smallest integer satisfying £ m = 1; we know that £ generates Q m . As 
£ G Q[e 2i7r/n ] too, Q m C Q». This does not preclude m > n. We need the 
following 

Lemma 2 77ie multiplicative group of elements of finite order in Q„ is cyc/icJl 

This is because given two elements £, £' with orders m, ml it is possible to con- 
struct an element of order lcm(m, ml). In other words, the roots of unity in Q„ 
have a maximum order, which is the 1cm of all possible orders. 

Let us call again m this maximal value, to prove Lemma 1 we need to prove 
that m — n or In. Now, any element £ of Q n which is a root of unity must satisfy 
f m = 1. 

This is true in particular when £ is the primitive n th root e 2 ' 71 "/™; hence m is a 
multiple of n, which entails Q„ C Q m - Finally Q„ = Q m . 

n I m and 



Letting <p(n) = dim(Q n /Q) stand for Euler's totient function, 



tp(n) = cp(m) 



( 1 \ I n for n even 

As <p{n) — n I 1 J , the only possibility is that m = < 

p| n;p prime V V J I 2n for 71 odd 

This proves the Lemma[T] and hence the Proposition: all eigenvalues of U are n 

71 + 1 

or 2n th roots of unity. Let us clarify the case of odd n: e iw/n = -(e 2i7r/n ) 2 
and we do have Q n = Q2 n - So we can rephrase the Lemma: in the odd case, 

C = ±1. 

Remark 1 At this point, U could be constructed as a polynomial in the elementary 
circulating matrix J (as all other circulating matrices) U = P{J), where P is the 
interpolating polynomial that sends the Fourier coefficients of J, i.e. the e 2lk7T / n t 
to the Fourier coefficients chosen for U. Such a construction is easy, with the 
Lagrange polynomials associated with the e 2lfc7r /™ ) as p i s a linear combination 
of these polynomials with coefficients that are precisely the Fourier coefficients of 
the desired u. 



For instance Q3 = ■ 

3 It is perhaps not obvious that this group is finite, and indeed the group of elements of Q n with length 
one is not; this holds because for large m the dimension of <Q m /Q> exceeds that of Q n /Q, equal to ip(n), 
cf. below. 
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3.2 Proof of Thm.|2] 

The possibilities of mapping the n th roots of 1 to m th roots of 1 can be somewhat 
reduced by noticing that U is a rational polynomial in J, and such a polynomial is 
stable under all field automorphisms of Q„ if we use the following characterization 
from Galois theory: 

Lemma 3 Any object (number, vector, polynomial, matrix) with coefficients in Q n 
is rational-valued iff it is invariant under all Galois automorphisms of the cyclo- 
tomic extension Q n over Q. 

We mention without proof either the structure of the Galois group. These two 
results can be found in any textbook on Galois theory. 

Lemma 4 Any field automorphism of the cyclotomic extension Q„ over Q is de- 
fined by $ fe (e 2l7r/n ) = e 21 ^ 1 " for some k £ Z*, the group ofinvertible elements 
of the ring Z n , e.g. for any integer k coprime with n. 

This is enough to define $t (x) for any x £ Q n , as x — J2aje 2ijn/n with rational 
a,j 's, and hence &k(x) = ^2 aje 1 " 1 ^ ^ n . For instance when n = 12, the four dif- 
ferent automorphisms $^ are defined by the possible images of e 2l "V 12 — e l7T ^ s , 
namely e iW6 ,fc € {1,5,7,11}. Their group is isomorphic with the multiplica- 
tive group Z* r2 = {1, 5, 7, 11}. 

If is such an automorphism, notice that = £, k for any n th root £ of 

unity (with one exception: 1) = — 1 Vfc £ Z£). If n is odd and £ is a 2n th 

root but not a n th , then — £ is a n th root, and hence 

£k for odd k's 
— (,k for even k's 

For instance $2 (£) = — £ 2 for such £, 

So from Lemma[3] we will have that U G M n (Q) if U is invariant under all 
the k G Z*. 

Now we can prove Thm[2] beginning with the case of n even. 

Consider the eigenvector X 3 = (1, e 2ijn/n , e 2i2 ^/™ 5 . . . e 2y(n-iW»y for 
the eigenvalue £j of f/ (for matrix J, the eigenvalue is of course e 2 ^/™). The 
prime signals that we consider Xj as a column. We have &k.(Xj) = Xjk by 
direct computation. The case j = is straightforward, as the eigenvector is real 
valued, so must be the eigenvalue, i.e. £0 = ±1. We exclude this case from now 
on. 

We assume that = U (i.e. that f7 is rational valued). Applying the field 

automorphism $fe to the equation 

UXj = ZjXt yields <k k (U)* h {Xi) = UX k] = = 9 k (^)9 k (Xj) = 

4 The coefficients of P can be read on the first column of U. 



MO = -(-£)* 
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Hence 

= # = (to 

for all j / and all k £ Z^. 

Now for the reciprocal. Assume the above equation (jj) between the eigenval- 
ues. We choose one Galois automorphism, $fe (for some k coprime with n). Let 
us apply <E>fc (U) to any eigenvector Xj ofU: notice thatX, = $k(X k -ij) where 
k~ 1 j is computed modulo n. Hence 

$ k (u)x 3 = * h (ux h -i s ) = Skfa-^Xt-ij) = ML-^)Mx k -i,) 

= £ k -ijX kk -ij because $fe raises any root of 1 to the k th power 

= (,jXj by our assumption on the eigenvalues 

We have proved that 3>fc(?7) does the same thing as U on any eigenvector; these 
make up a basis, hence = U, i.e. U is rational valued. 

It remains to discuss the case of odd n. We still get the equation $fc(£j) = 
if U is assumed to be invariant under $fc . 

If £ is a n th root of unity, the computation is identical. 

If i 2n = 1 but C / 1, then (-£)" = 1 and hence = -$*(-£) = 

— (— = — £ fc for even fc and = £ fc f° r °dd k. The computation above 

still yields = 4 > fc(Ci) = f J f° r odd fc, and we have also the new case £j k = 
— for even k. 

Say fe = 2, and £i = £ with f? n = 1 ^ then = -£ 2 ,f 4 = 
— £ 4 , • • • £2™ = — £ 2 • 2 has a finite order in Z* , hence for some m, £2™ = £1. 
We get an orbit of m eigenvalues which are all 2n th roots of unity, e.g. O = 

{fl,£2,f4,£8 ■ ■ ■ }■ 

Say now that k = 2 v k' ,k' odd and coprime with n. The formula (jj) is then 
valid and yields = £|« ■ So is determined when O is known. Notice that 
£fc will never be a n th root (because 2 and k' are coprime with n): either all the 
eigenvalues [with even index] are n th roots, or none (except of course £0 = ±1). 

The reciprocal is similar to the even case: it is identical when the eigenvalues 
are of order n (at most); and if £1 has order 2n then the values of £fc that we have 
obtained enable to satisfy the relations $k(U)Xj = £jXj for all j, k k coprime 
with n), so that <&/t(f7) is identical to U i.e. U is rational- valued. This ends the 
proof. 

3.3 Proof of Thm.|3] 

We make us of the conditions found in Thm. [2] 

The whole set of eigenvalues is thus determined if we know £j for a subset of 
representatives j of all orbits under multiplication by elements of Z* (so called 
associated elements in the ring Z n ). We can specify the smallest representatives: 
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Lemma 5 Any element j £ 7L n is associated with a divisor of n, i.e. 3 A: £ 

= gcd(n, j). 

(We identify integers and classes modulo n when the distinction is irrelevant). 

This stems from the Bezout identity (in Z): for some k, £, kj+£n = gcd(n, j). 
After division by gcd(n, j) we see that k and n are coprime. But modulo n, kj — 
gcd(n, j), qed|f| 

So it is sufficient to specify £j when j is any divisor of n. We will need a last 
Lemma, which seems interesting in its own right@ 

Lemma 6 The set of differences A n = Z* — Z* = {a — b, (a, b) £ (ZJj) 2 } is 
Z n when n is odd, 2Z n when n is even. 

It is straightforward for n prime, for n an odd prime power, and we notice that 
when n = 2 m then 7L* n = odd numbers, so that A n — even numbers. The general 
case now stems from the Chinese remainder theorem. 

We now procede to prove the theorem. Remember that £o = ±1. 
3.3.1 When n is even 

In this case all eigenvalues are n th roots of unity. Let j be any strict divisor of n. 

• When n/j is even, we can produce k, k! £ Z* with k! — k = — £ 2Z from 

3 

Lemma|6] Hence (noting that k = k' mod n) 



which proves that ^ 3 = 1, i.e. is a power of e 



2ijiv j r 



• If n/j is odd (meaning that j contains the same power of 2 as n), then 

Lemma[6]only provides k' — k = — , and the calculation yields £^ = 1, 
..3 

i.e. is a power of e 1 ^/™, which ends the even case of the theorem. 



3.3.2 When n is odd 

The case when ^ is a n th root is identical to the n even (first) case, as from the last 

Lemma[6] we can again produce two elements k, k' £ V n such that k' — k — —. 

j 

and t;* _fc = 1 = Cj^ ■ So the spectral unit is determined when we have chosen a 
n/j th root of unity £j for each divisor j of n. 



5 For instance for n = 15 we have the orbits of equivalent elements 

(0), (1, 2, 4, 7, 8, 11, 13, 14), (3, 6, 9, 12), (5, 10) 

indexed by the divisors 1,3,5 and of course 0. 

6 Though 
the literature 



6 Though elementary, the result was unknown to the author and does not appear to be readily available in 
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Now assume that there is an eigenvalue £j which is not a n th root. Then — 
is a n th root, and (as n/j is odd) a similar calculation yields for k! — k = n/j, 
with fc, fc' € Z* , 

= = -ijk = (-o) fe = (-0) fc ' _ * 

hence — £j is again a n" 1 root of unity. 

This ends the proof of the odd case of Thm. [3] 
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